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Mass Exchange Processes with Input 


P. L. Krapivsky 1 
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We investigate a system of interacting clusters evolving through mass exchange and supplemented 
by input of small clusters. Three possibilities depending on the rate of exchange generically occur 
when input is homogeneous: continuous growth, gelation, and instantaneous gelation. We mostly 
study the growth regime using scaling methods. An exchange process with reaction rates equal 
to the product of reactant masses admits an exact solution which allows us to justify the validity 
of scaling approaches in this special case. We also investigate exchange processes with a localized 
input. We show that if the diffusion coefficients are mass-independent, the cluster mass distribution 
becomes stationary and develops an algebraic tail far away from the source. 


PACS numbers: 82.20.-w, 05.40.-a, 89.75.Fb, 68.43.Jk 

I. INTRODUCTION 

Exchange processes underlie numerous phenomena. In 
natural sciences, exchange plays an important role in 
droplet growth via evaporation and re-condensation [T], 
island growth in deposition processes [2], and phase or¬ 
dering EHS1- Exchange processes have been also applied 
to modeling segregation of heterogeneous populations [B], 
studying the distribution of wealth B, mimicking growth 
of urban populations [8j and aggregation behaviors in job 
markets !0 . 

The basic ingredient of exchange processes are clusters. 
We suppose that clusters are composed of a certain in¬ 
teger number of elemental building blocks (‘monomers’); 
clusters are thus labelled by their masses, other charac¬ 
teristics (e.g., their shape) are not explicitly taken into 
account. We also assume that in each exchange event, 
a monomer is transferred from one cluster to another. 
An exchange event can be separated into detachment 
of a monomer, subsequent transport to another cluster 
and attachment to it, so the comprehensive description 
of exchange processes can be complicated. We focus on 
the detachment-controlled limit when the time scale for 
transport and subsequent attachment is negligible. 

Denote by Aj a cluster of mass j, that is, a cluster 
which is made of j monomers. Symbolically, the mass 
exchange process can be represented as 

Ai + Aj — A Ai ±i + Aj^zi 

We assume that a cluster is equally likely to gain or to 
lose mass in an interaction. This implies that the reaction 
rate matrix (‘kernel’ in short) is symmetric: j = Kj :i . 

A cluster disappears when its mass reaches zero. Thus 
the number of clusters decreases by 1 with probabil¬ 
ity 1/2 when strictly one of the reacting clusters is a 
monomer; if both reacting clusters are monomers, one of 
the monomers disappears and the other becomes a dimer. 
Exchange processes with rather general classes of homo¬ 
geneous symmetric kernels, e.g., with kernels of the form 
Kij = i a j b + i b j a , have been studied in Refs. nnna. 

The density c m {t) of clusters containing m monomers 


at time t evolves according to the rate equation 

—ffj— = ^ ^ CjCj Kjj \6 m ,i+l T 6m j— 1 2(5 m y] (1) 

This equation assumes perfect mixing, or equivalently, 
absence of spatial correlations. Furthermore, the total 
number of clusters is assumed to be infinite; otherwise, 
even if we ignore the locations of clusters, the system 
will eventually condense into a single cluster and close 
to the condensation time fluctuations in number of clus¬ 
ters become important, so the deterministic rate equation 
framework looses its validity. Even for infinite systems in 
homogeneous settings, the rate equation framework be¬ 
comes inapplicable in low dimensions, particularly in one 
dimension (see [T21 H3] and references therein). 

If the reaction kernel is the product of reactant masses, 
Ki j = ij, the cluster densities c m (t) satisfy 

dc 

—AA = M[{m + l)c m+ i - 2 mc m + (m - l)c m _i] (2) 

and since the mass density M = '}2j >1 jcj remains con¬ 
stant due to mass conservation, Eqs. H are linear. The 
linearity of rate equations distinguishes the product ker¬ 
nel from other kernels. If the system is initially composed 
of monomers, c m {t = 0) = <5 m> i, the solution reads 

= (i + *)">+! (3) 

The corresponding total cluster density is 

AT = ^c i = (l + t)- 1 (4) 

i> i 

The cluster distribution ([3]) is the simplest exact result in 
the context of exchange processes BE]; it was originally 
discovered in earlier studies of birth-death processes (see 
m and references therein). A few other exact solutions 
for exchange processes with more complicated kernels can 
be found in BEJ E3II6]. 

Exchange processes help in classifying and interpreting 
emergent phenomena in various complex systems W- 
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Pure exchange process may provide an idealized descrip¬ 
tion of binary interactions underlying the dynamics of 
such complex systems. In addition, complex systems are 
usually open and driven. Perhaps the most ubiquitous 
example of an external drive is input; e.g., it is a key 
element in deposition processes and it naturally arises in 
phase ordering. The goal of this paper is to study the 
interplay between exchange and input. 

Pure exchange processes resemble aggregation pro¬ 
cesses m- In the framework of aggregation, the influ¬ 
ence of input has been studied. The outcome was some¬ 
what surprising: For basic reaction kernels, aggregation 
processes with homogeneous input are simpler than cor¬ 
responding pure aggregation processes, viz. in the most 
interesting large time limit the dominant part of the clus¬ 
ter mass distribution acquires a stationary form (see e.g. 
[T5l;16l (181121] 1. There are exceptional kernels and when¬ 
ever the stationary form is never reached, the emergent 
behaviors tend to be subtle (often they are character¬ 
ized by more than one mass scale [2T1 122] ). We shall 
show that in contrast to aggregation, exchange processes 
with spatially homogeneous input do not reach steady 
states. The emergent evolutionary behaviors, however, 
are usually qualitatively the same as for exchange pro¬ 
cesses without input mi. 

Our second goal is to study exchange processes with a 
localized source. A classification of possible behaviors in 
such processes appears unattainable—the behaviors de¬ 
pend on the reaction kernel, on the (mass-dependent) dif¬ 
fusion coefficients, and mathematically one has to study 
an infinite system of coupled non-linear partial differen¬ 
tial equations. In the idealized situation when the dif¬ 
fusion coefficients are mass-independent, however, the 
models with a localized source are partly tractable. More 
precisely, the cluster mass distribution becomes station¬ 
ary, at least in the most relevant three-dimensional set¬ 
ting, and this feature greatly simplifies analysis. 

Our analysis of mass exchange processes with input 
relies on scaling techniques. The emergence of scaling is 
difficult to prove, so it is useful to have exact results at 
least for simplest kernels which can be compared with 
outcomes of scaling analyses. It turned out that the 
model with product kernel remains solvable when we add 
homogeneous input. A few basic densities are still given 
by neat formulas, especially in the most natural setting 
when (i) only monomers are injected, and (ii) the system 
is initially empty. We will show that in this situation the 
total density reads 


1 ln Vt 2 + 2+ t 
Vt 2 + 2 n Vt 2 + 2 - t 


while the density of monomers is 


(5) 


ci = 


i n vgT2 + t (6) 


t 2 + 2 ( t 2 + 2)3/2 JpT2 - t 


the exchange process with product kernel and homoge¬ 
neous input of monomers. This is a unique kernel for 
which an exact approach (based on the generating func¬ 
tion technique) has been successful so far; another exactly 
soluble kernel (Sec. IV C) is degenerate as it only cou¬ 
ples monomers and dimers. We are mostly interested in 
long time behaviors, however. Such behaviors can be ex¬ 
plored using scaling techniques 

unna is]. Additionally, 
scaling techniques are applicable not just to the solvable 
product kernel, but to numerous kernels, even families of 
kernels. In Sec. |III| we illustrate scaling techniques for 
exchange processes with three basic kernels: (i) product 
kernel Ki j = ij, (ii) constant kernel Kij = 1, and (iii) 
sum kernel Kij = i + j. In Sec. IV we analyze exchange 
processes with a one-parameter family of kernels, namely 
the generalized product kernels K t j = ( ij) x . We show 
that scaling methods are applicable to the generalized 
product kernels with homogeneity index A < |. In the 
| < A < 2, an infinite cluster (‘gel’) is formed at a finite 
time, while for A > 2 there is instantaneous gelation. 
In Sec. [TV] we also briefly discuss the extreme values of 
the homogeneity index: A = — oo which is exactly sol¬ 
uble and A = 2 showing a multi-scaling behavior. In 
Sec. [V] we study exchange processes driven by a localized 
source. We show that when the diffusion coefficients are 
mass-independent, the densities attain a steady state, at 
least in the most relevant three-dimensional setting, and 
we determine these densities for a few basic kernels. We 
also probe models with diffusion coefficients decaying as 
(mass) -1 and show that a stationary scaling mas s distri¬ 
bution emerges when A < We close (Sec. VII with a 
short discussion. 


II. EXACT ANALYSIS 

Here we investigate the mass exchange process with 
product kernel and homogenous input using exact tech¬ 
niques. We shall assume that only monomers are in¬ 
jected. This is the most natural choice: The same ele¬ 
mental mass which is transferred in each exchange event 
is also added. The details of the input, however, do not 
affect qualitative behaviors (as long as clusters of only 
small masses are added). We shall also assume that the 
system is initially empty, 

emit = 0) = 0, (7) 

so that the mass density is equal to time, 

J2i c j( t ) = t , ( 8 ) 

i> i 

if only monomers are injected and if the monomer flux is 
set to unity. 

In the case of the product kernel, = ij , we thus 
arrive at the infinite set of rate equations 


The remainder of this paper is organized as follows. 
In Sec. [I] we derive © © and other exact results for 


dCm 

dt 


= <5r7 


+ t[(m+l)c m +i -2rnc m + (ro-l)c m _i] (9) 
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Equations |9j) are linear. In this respect, they are sim¬ 
ilar to Eqs. (pi describing the evolution without input. 
In contrast to Eqs. (pj), however, we now have an in¬ 
homogeneous system of equations with time-dependent 
coefficients. It is still possible to solve (pj) subject to the 
initial condition (pj). Multiplying (p| by z m and summing 
over all to > 1 we reduce an infinite set of rate equations 
© to a single partial-differential equation 


1 d_ 
t dt 


{l-zf 


d_ 

dz 


e 


t 


- ci(t) 


( 10 ) 


for the generating function 



C(M) = 




m> 1 


Transforming variables ( z , t) to variables 

1 


£=-- 
? 2 1-z 


i 2 1 

^=2 + l^z 


we recast (10) into 

oe 


1 - 


2 — = 




dt, Vt + V 


- Cl 


(Vt + v) 


(ii) 


( 12 ) 


(13) 


Integrating (13) yields 


<'ie,n) = J 


ft du 

'1 

1 

J-n V u + V 

2 

T] — U 


f 


du Cl (y/f. 


(14) 


The choice of the lower integration limit in ( |14[ ) ensures 
6|| = _,j = 0, which is equivalent to G(z, t = 0) = 0, i.e., 
to our choice of the initial condition, Eq. (Hi- 

Computing the first integral in( |14| ) and returning to 
the original ( z , t) variables we obtain 


C (z,t) =t— dvvc\{v) 

Jo 

y/t 2 2^ + t 


In 


y/t 2 + 2C y/t 2 + 2( -t 


(15) 


where ( = (1 — z)~ 


According to the definition (|ll|), the generating func- 

1 

1 


tion must vanish for z = 0. Therefore 

rt 


/ dv vci(v) =t — 
Jo 


Vt 2 


l n vg+2 + t (16) 


Vt 2 


2 — t 


leading to the announced expression p]) for the monomer 
density (see also Fig.p]). 


Using (161 we can re-write (15) as 
1 


G(z,t) = 


In 


Vt 2 + 2+ t 


Vt 2 + 2 y/t 2 + 2-t 
1 \/t 2 + 2C + t 

y/t 2 + 2C H yjt 2 + 2( - t 


(17) 


FIG. 1: Evolution of the total density N(t) and the monomer 
density ci ( t ) in the mass exchange process with product kernel 
and homogeneous input. The top curve shows the total den¬ 
sity given by Eq. |5|; the bottom curve shows the monomer 
density given by Eq. §. 


Since 6 (z = 1 ,t) = N(t) we obtain the total cluster den¬ 
sity by specializing (15) to z = 1. This gives pj shown 
on Fig. p] 

The generating function <0 encapsulates all cluster 
densities. Expanding C (z,t) in power series in 2 one can 
extract explicit expressions for cluster densities. General 
expressions for cluster densities are similar to © , viz. 


tPm(t ) 


Qm(t) 


ln V tH + t (18) 


(t 2 + 2) m (t 2 + 2) m + 1 / 2 Ut^+2 - t 


where P m (t ) and Q m (t) are even polynomials of degree 
2 (to — 1). For instance, 


Pi = 1 1 1 

*-!*• 


1 

2 ’ 


5,2 

6 l 


1 

2 ’ 


p — 1,6 _ 13,4 _ 31 u 
-T4 — 4 t 12 l 9d i 


24 


Qi — t 2 + \ 

Qo = t 4 +1 2 + \ 

|, Q 4 = t 6 + lt 4 + %t 2 + l 


III. SCALING ANALYSIS OF THREE BASIC 
KERNELS 

In this section we employ scaling techniques. We start 
with the product kernel for which one can justify scaling 
using an exact solution from Sec. [XT] We then turn to two 
other simplest kernels, the constant kernel and the sum 
kernel, for which scaling appears to be the only available 
analytical technique. 


A. Product Kernel 

Here we re-consider the product kernel, Kj j = ij, but 
use scaling instead of exact methods. The typical mass is 
expected to grow indefinitely. For large masses, we can 
treat m as a continuous variable and approximate the 
system of rate equations © by the partial differential 
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equation for the density c = c m (t): 


dc d 2 (mc) 

dt dm 2 


(19) 


Balancing both sides of this equation yields \ ~ —, so 
the typical mass grows according to to ~ t 2 . 

To establish the asymptotic behavior of the density of 
monomers we sum all Eqs. § and obtain the exact rate 
equation for the total cluster density: 


dN 

jT = 1 - tcl 


( 20 ) 


This equation suggests that 


lim fci = 1 

t—> OO 


( 21 ) 


The rate equation for the density of monomers 

dc i „ . , 

— = 1 - 2f[c, - c 2 ] 


( 22 ) 


suggests that lim^oo t(c± — C 2 ) = 1/2, which in conjunc¬ 
tion with (21) give 


lim tco = - 

t—> 00 2 


(23) 


Similarly analyzing the rate equation for the density of 
dimers we get linp^^ tc% = |, and generally for any fixed 
mass k we have 


lim tc k = y 

t-¥ 00 ft 


(24) 


Using (24) together with the growth law, (to) ~ t 2 , of 
the average mass we arrive at the following scaling form 
of the cluster distribution: 


c m (t) = - 
tm 


M = 


TO 

i 2 


(25) 


More precisely, the scaling form (25) is expected to hold 
in the scaling limit 


t 


u = — = fixed 
^ t 2 


(26) 


The small mass behavior (241 and the mass conservation 
law ([s]) imply the following properties of the scaled mass 
distribution: 


F( 0 ) = 1 

/»oo 

/ dnF(n) = 1 

J 0 


(27a) 

(27b) 


The scaling form ([25]) together with (27a) lead to a rather 
unusual asymptotic behavior of the total cluster density: 


N ~ - In t 
t 


(28) 


F 



FIG. 2: The scaled mass distribution F(fi) given by Eq. (31) 
arising in the mass exchange process with product kernel and 
homogeneous input. 


To determine the scaled mass distribution we plug the 
scaling form (25) into the governing partial differential 


equation (191 and find that the scaled mass distribution 


satisfies an ordinary differential equation 


F" + 2F' + - F = 0 


(29) 


where the differentiation with respect to /r is denoted by 
prime. Seeking solution in the form -F(/.i) = e -A1 <£’(/i) we 
eliminate the term with first derivative: 


4)" + - 1^ $ = 0 


(30) 


WKB theory shows that there are two possible large /i 
asymptotic behaviors, /u -1//2 e M and Only the 

latter asymptotic is acceptable. The boundary condition 
F{ 0) = 1, equivalently 4>(0) = 1, then leads to the unique 
solution, 4> = n[K 0 (fj.) + Ki(n)\, in terms of the modified 
Bessel functions K 0 and Ki. Thus 


F(/x) = ne M [A'o(m) + K i(/x)] (31) 


(see also Fig. [ 2 ]. 
write pi]) as 


Recalling that K\ = — A/j, we can re- 




(32) 


Plugging (32) into (27b) and perfor ming the integration 
by parts we see that the validity of (27b) is equivalent to 

(33) 


/»oo 

/ d/ie _M /\ 0 ( n) = l 

Jo 


This identity can be derived e.g. from the integral rep¬ 
resentation of the modified Bessel function Kq: 


K 0 (jm) 



du 


y/u 2 — 1 




(34) 


























5 


We also mention asymptotic behaviors of the scaled 
mass distribution. In the small mass limit (/u —> 0) 


(35) 


F(n) — 1 — /z ^hi — — 1 — 7 

and in the large mass limit (/z —> oo) 

~ sJ'l-Kne^ 


(36) 


For the product kernel we can use exact results to ver¬ 
ify scaling predictions. The exact prediction ([ 6 ]) for the 
monomer density gives 

Ci = — + — [2 In t — 2 + In 2 ] + ... (37) 


when t 1 , so the scaling prediction ( 21 ) agrees with 


the exact prediction in the leading order. Similarly, the 
exact prediction © or the total density, 

N ~ - [2 In t + In 2] + [— In t + 1 — In 2] + ... (38) 

for t 1 , has the leading order term which agrees with 
the scaling prediction ( 281 . 


To confirm the emergence of the mass scaling (26) and 


deduce the scaled mass distribution it is convenient to re¬ 
write the exact prediction ( fl7| ) for the generating function 
in the form 


5^(1 - z m )c m (t) = 


1 In V /f2 + 2 ^ + f (39) 


1>1 


yjt 2 + 2C yjt 2 + 2C -t 


The scaling regime emerges when £ ~ t 2 . Recalling that 
£ = (1 — ,z) _1 , we write 


2 = 1 “? 


(40) 


and take the f —> oo limit while keeping s finite. Com¬ 


bining (39) and (40) we see the consistency of (25)—(26). 


Furthermore, plugging (40) and (25) into (39) we obtain 
an integral equation 


r - (i - e ~nm = in 

Jo f 1 \ s + 2 


y/ s T 2 + y/s 

y/ s + 2 — yfs 


(41) 


for the scaled mass distribution. One can verify that the 
scaled mass distribution which we established using scal¬ 
ing methods is indeed the solution of the integral equa¬ 
tion d41). 


B. Constant Kernel 

For the constant kernel K, j = 1, the rate equations 
read 

dc 

—jj~ = l$m,l + N[c m + 1 — 2 C m + C m _ i] (42) 


Summing (42) we obtain 


dN „ Ar 

s = l - Nci 


(43) 


For large t and m we can approximate the system of 
rate equations (42) by the partial differential equation 


which happens to be the diffusion equation 

d 2 c 


dc 


dm 2 


where we have use an auxiliary time variable: 


= f dt'N(t') 

Jo 


(44) 


(45) 


The appropriate scaling form for the cluster mass dis¬ 
tribution is 


i ( t ) = 


= -r-l/4 




M = -p (46) 


The form of the scaled mass, fi = m/y/r , is obvious from 
(44). The time-dependent pre-factor in (46) was chosen 


to ensure the validity of the asymptotic, Nc\ —> 1 when 
r —► oo, which is implied by Eq. (43). Indeed, using (46) 
we find that the cluster density N{t) = Z) m > i c niJj) 
grows according to 


N = t 1 /^ / d/iF(jj,) 
Jo 


(47) 


Therefore iVci is asymptotically constant and it equals 
to unity if 


/»00 

F( 0 ) / dnF(») = 1 
Jo 


(48) 


Hence N = [F(0)] 


-WV4 


. Plugging this into (45) we get 


t 4F(Q) _3/4 


Using (|46| we compute the mass density 

/»oo 

55 JCj (t) = r 3/4 / d/j.nF(fj,) = t 
Jo 


3> 1 


Comparing ( |49[ ) and ( |50[ ) we obtain 

/»oo 

cfyz [iF(fj.) 


4F(0) 


(49) 


(50) 


(51) 


To complete the analysis we must find the scaled mass 
distribution F(/z). Plugging (46) into (44) leads to an 
ordinary differential equation 


F" + 4/zF 7 + 4 F = 0 


(52) 


Demanding that F(/z) decays sufficiently rapidly we ar¬ 
rive at a one-parameter family of solutions 

F = C y/Ji exp(—/ z 2 /8) A" 1/4 (/z 2 /8) 


(53) 
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involving the modified Bessel function K\i 4 . 


To fix the amplitude C we use relations (48) and (51). 


In computations, an integral representation of the mod¬ 
ified Bessel function K 1 i 4 proves useful. This integral 
representation is obtained by specializing 


K v {z) = 


r(- + !) 


IT©'. 


du 


y/'U 2 — 1 


(u 2 -ir e - 


to v = 1/4. Using above integral representation one can 
compute the two-parameter family of integrals, 


POO 

K(0,v):= dz z^~ l e~ z K u {z) : 

Jo 

with arbitrary /?, v satisfying /3 > v. One finds 


(54) 


K(P,v) = 


iP + v) 


du 


(u- 1)^-3 


2-r(z/ + i) A (« +1)^+1 

0rT(/3 + z/)r(/3 - z/) 


2/3 r (/3 


(55) 


As a check we note that if (1,0) reduces to the integral 
in (33) and the right-hand side of (55) is equal to unity 

for [Pi v) = (1,0). 

Using (53) and writing /i = y/8z we express the in¬ 
tegrals appearing in (48) and © through the integrals 
K(j3,v), namely 


POO 

/ d^F^) = C-2 5 / 4 K(l\) 

Jo 

POO 

/ d / i / iF( M ) = C'-2 11 /4A'(|,|) 
Jo 


(56) 


Combining (56) with (551 and adding F (0) we get 

F( 0 ) =C W ^ 


r (!) 

L d >‘ F M = c ¥§j 

POO I7y 

/ dnnF(n) = C-j^~ 

Jo 1 ( 4 ) 


4 

X 3 


(57a) 

(57b) 

(57c) 


Equations (57a) and (57c) show that fl5lfo is manifestly 


satisfied. Using (57a) and (57b) we see that to ensure 


the validity of (48 1 we must choose 

1 


C = 


2 5/4, 


= 0.237212499... 


(58) 


The asymptotic behaviors of the scaled mass distribu¬ 
tion are 


f(h) - m = 


2 3 / 4 r (|) 


= 0.860039988... 


when (j, 0 and 


F(/x)~ 


2 1 / 4 v 7 r ' 


0 -M 2 /4 


(59) 


(60) 


when /r —> 00 . The latter formula implies that 


i( T ) - 


1 


2 1 / 4 y / m 


exp 


m 
4 T 


(61) 


for m yfr. 

Finally, we give asymptotic formulas for the total den¬ 
sity and the monomer density through the original time 
variable: 


n = At 1/3 , Cl = a-H~ 1/3 

where the amplitude is given by 


(62) 


A = 


6 1/3 [r(f)] 

7T 2 /3 


4/3 


= 1.1108674179... 


C. Sum Kernel 

For the sum kernel K t j = i + j, the rate equations 
read 


dCm 

dt 


= N[(m + l)c m+ i - 2 mc m + (to - l)c m _i] 


T 2 C m C m _i] + ^m,l 


The total cluster density satisfies 
dN 

-dT = 1 - tCl 


(63) 


(64) 


leading to the same asymptotic relation ( 21 ) for the 


monomer density as in the case of the product kernel. 

In the scaling regime, the cluster density c = c m (t) 
obeys the partial differential equation 


dc d 2 (mc) d 2 c 

dt dm 2 dm 2 


(65) 


Balancing the left-hand side with the first term on the 
right-hand side gives m ~ Nt, while balancing with the 
second term leads to to ~ t. These estimates for the 
typical mass are comparable, and thence both terms on 
the right-hand side remain asymptotically relevant, if N 
saturates at a certain finite value N a0 . Below we compute 
this value and find 


lim N (t) = iVoo = 1 

t—> OO 


( 66 ) 


but now we merely suppose that the cluster density 
reaches a finite value N^. Since to ~ t, the cluster mass 
distribution is sought in the scaling form 

1 TO 

Cm{t) = -F(n), 


Combining (65) and (671 we arrive at 


F + nF' + [F + NocixF]" = 0 


( 67 ) 
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Integrating once we obtain /xF + [F + -ZVoo/xF]' = 0, or 
equivalently 


[1 + iV 00 /x]F , + [IVoo + /x]F — 0 


where r is an auxiliary time variable 


r = 


dt' 


( 75 ) 


Integrating this equation gives a family of solutions 


F = 


C 

(l + AW) 1 - 1 /^ 


g-M/Wc 


( 68 ) 


with two un-determined constants: C and -/V^. 

We should impose three conditions on the scaled mass 
distribution: 


F( 0 ) = 1 

/»oo 

/ dfi F(fi) = 

Jo 

/»oo 

/ d/x/xF(/x) = l 
Jo 


(69a) 

(69b) 

(69c) 


Equation (69a) ensures the validity of (21); Eq. (69b) 


follows from lim^oo X)m>i c m(t) = Noo ! the last condi¬ 
tion, Eq. (69c), is the conservation law X) m >i c m(t) = t. 


Equation (69a) obviously leads to C = 1. Further, one 


can verify that the only choice of N x which agrees with 
Eqs. (|69b|) (|69c|) is N, x = 1 as it was stated in (66|). Thus 


(68) simplifies to 


F = e -/i (70) 

IV. GENERALIZED PRODUCT KERNELS 

In this section we consider exchange processes with 
generalized product kernels, K li0 = ( ij ) x , and homoge¬ 
neous input of monomers. The rate equations read 

dc 

-dtt = M A [(TO+l) A c m+ i-2m A c m + (rn-l) A c m _i] (71) 

for 77i > 1, while the monomer density satisfies 

^ = l + M A [2 A c 2 -2 Cl ] (72) 

Here M\ is the A th moment of the cluster mass distribu¬ 
tion: 

M x (t) = ^2 mXc m(t ) (73) 

m> 1 


Note that the cluster density is the zeroth moment and 
the mass density is first moment of the cluster mass dis¬ 
tribution: N = Mq and M = M\. 

For large t and m we replace the system of rate equa¬ 
tions ( |7l| ) by the partial differential equation 

dc = d 2 (m x c) 

3t dm 2 1 ’ 


We seek the cluster mass distribution in the scaling 
form 


C m (r) = r F(n), d = — 


(76) 


The form of the governing equation (74) fixes the expo¬ 


nent a which determines the growth of the typical mass: 

1 


a = 


2-A 


(77) 


To find another exponent let us establish the growth law 
of M\. Plugging (76) into (73) we get M\ ~ T (i+x)a-b 


Combining this result with (75) yields 


t - 


1+b— (1+A)a 


(78) 


Computing the mass density gives another relation be¬ 
tween t and the modified time r: 


t = M(t) ~ r 


2a—b 


Relations (78) and (79) are compatible only when 


b = 


1 + 2A 
2(2 - A) 


(79) 


(80) 


The expressions ([77f and ( |80| ) for the scaling expo¬ 
nents already tell us that the above analysis is inapplica¬ 
ble when A > 2. For A > 2 the mathematical framework 
becomes ill-defined due to instantaneous gelation pTj . 

Instantaneous gelation is a rather counter-intuitive 
phenomenon implying that an infinite cluster (‘gel’) nu¬ 
cleates at time t = +0. Instantaneous gelation was orig¬ 
inally discovered in the context of aggregation [23H2S], 
see for recent work. Instantaneous gelation also 

arises in addition processes [213 where it was rigorously 
proved [30]. Instantaneous gelation was mostly studied in 
the mean-field framework, but it may also occur in one 
dimension HUSH where spatial effects are crucial and 
the mean-field rate equation approach is inapplicable. 

In the range | < A < 2, one can employ scaling anal¬ 
ysis using the auxiliary time and considering the usual 
r —> oo regime. One cannot use the original time vari¬ 
able, however. This can be appreciated by using (77), 
(80) to re-write (79) as £ ~ r ( 3 - 2A )/( 4 - 2A ) indicating that 


consistent behaviors occur only when A < |. In the 
| < A < 2 range the system actually exhibits gelation, 
namely the gel forms in a finite time. Interestingly, gela¬ 
tion is complete meaning that the entire mass is sud¬ 
denly transformed into gel. The analysis of the behavior 
in the gel regime is essentially the same as in un-driven 
exchange processes with the same kernels m3- Indeed, 
the most interesting behavior arises close to the gelation 
point, and the mass density is finite as in the un-driven 







































case. Therefore let us focus on the A < § range where 
clusters grow indefinitely and the emergent scaling be¬ 
havior differs from the un-driven case. 


Instead of the scaling form (76 1 , let us use an equiva¬ 
lent, but slightly more convenient, scaling form 

Cm(T) = T- a/ 2 m- x ${n), /x=^, a= J~x ( 81 ) 

Plugging this scaling form into the governing equation 
(74) we arrive at the ordinary differential equation 

+a + !$] = 0 (82) 

Demanding that >1>(/t) decays sufficiently rapidly we ar¬ 
rive at a one-parameter family of solutions 

$(g) = CiFi(a/2,l^a;-aV / “) (83) 

in terms of the (confluent) hypergeometric function 

T(/3) ^ T(a + n) z n 


iFi(a,/3;z) = 


E 


r (a) ^r(/3 + n) n! 


Using (El) and (|83|) we compute the A th moment: 


M x = CT a/2 h 


with 


L = a 


pOO 

/ dzz 0 '- 1 iFi(a/ 2,1 — a; — a 2 z) 

Jo 


(84) 


(85) 


Similarly we compute the first moment, the mass density: 

M = t = CT 1 ~ a/2 I 2 ( 86 ) 

with 


h = a 


p OO 

/ dz iFi(a/ 2,1 — a; — a 2 z) 

Jo 


(87) 


bined with (84) to give 


Equation (75) tells us that 44- = 1 /M\ which is com- 


= T -a /2 _L 


dr 


CL 


Differentiating ( 86 ) with respect to r we obtain another 
expression 

^ = T -“/ 2 (1 - a/2)CI 2 
dr 

for the same derivative. These two formulas are compat¬ 
ible when 

1 


C = 


\/(l — a/ 2 )Ji / 2 
Note also the behavior of the total cluster density 


( 88 ) 


N = C7n t (!-2A)/(4-2A) = (jj 


CL 


1-2X 

3-2X 


with 


In = a 


r°° dz 

/ -iF, (a/ 2,1 -cr,~a 2 z) 
Jo z 


(89) 


(90) 


Let us look more carefully at a few particular values of 
the homogeneity index. 


A. A = | 


The scaled mass distribution substantially simplifies 
when the homogeneity index is A = In this case the 
parameters of the confluent hypergeometric function co¬ 
incide (| = 1 — a = 3 ), so the function becomes expo¬ 
nential. More precisely, 


$(/x) = Cex p[-(f)V/ 2 


(91) 


with 


' 9 \ 1/6 r 

c = 'j) 




1-1/2 


= 0.8031990959... (92) 


The total cluster density in this case saturates at 

1 3/2 


N °°~ 3 


r (f) 


= 1.87970144... 


(93) 


B. Exchange process in d dimensions 

Consider an exchange process supplemented by input 
of monomers in d spatial dimensions. We assume that 
clusters are spherical (say because of the surface tension) 
and that in an exchange event a monomer on the surface 
of one cluster migrates to the surface of another clus¬ 
ter. This is similar to the model proposed by Schelling 
[ 6 ] to account for segregation. Another interpretation 
is an Ising model with infinite-range zero-temperature 
Kawasaki spin-exchange dynamics m- The precise re¬ 
lation between all these models is rather subtle, e.g., 
there is no clear surface tension in the simplest versions 
of the Schelling and Ising-Kawasaki models, the shape 
of clusters apparently reflects the lattice structure El- 
Coming back to our exchange process we note that the 
corresponding kernel is the product of the surface areas 
of the clusters participating in exchange: = cqoj. 

Expressing the area of the cluster and its mass through 
the radius, a ~ R d ~ l and m ~ R d , we conclude that 
<7 m ~ to 1 ” 1 /^, so this exchange process is described by 
the generalized product kernel with A = 1 — 

In one dimension A = 0, so the exchange process has 
constant kernel (Sec. |IIIB |. In particular, the scaled mass 
distribution is given by (53) and (58), the total cluster 
density grows as t 1 / 3 , Eq. (62). In two dimensions A = 
2 , so we recover the exchange process described in the 
previous subsection |IV A| The scaled mass distribution in 
two dimensions is given by (91) (92) and the total cluster 

In three 


density saturates at a constant value, Eq. 
dimensions A = |, the scaled mass density becomes 

$(/x) = C 1 F 1 (f,i;-^ 4 / 3 ) 


(94) 


and the total cluster density decays as t 1 / 5 . Generally 
in d dimensions: N ~ /( 2 ~ d )/( 2 + d ). 
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C. A = -oo 

For this smallest value of the homogeneity index, the 
scaling approach is no longer applicable. Yet the analysis 
of this case is simple since (i) only monomers participate 
in exchange, and (ii) in addition to monomers injected 
into the system, only dimers can be formed. 

The rate equations give a single non-trivial equa¬ 
tion 


dc 2 

dt 


= Ci 


(95) 


while (72) reduces to 


dc\ 

dt 


= 1 - 2 cf 


(96) 


Solving (95) (96) subject to ci(0) = C 2 ( 0 ) = 0 we obtain 

Ci = —j= tanhT, C 2 = — -j= [T — tanhT] (97) 

V2 2y/2 


with T = ty/ 2. 


D. A = 2 


The value A = 2 of the homogeneity index is another 
extreme: It separates the | < A < 2 region where the 
gel is formed from the A > 2 region where instantaneous 
gelation occurs. The simplest way to convince oneself in 
the emergence of a gel, and to establish the gelation time, 
is to look at the moments of the cluster mass distribution. 
In the particular case of A = 2 the equations for (integer) 
moments are recurrent , and hence solvable. For instance, 


d A/2 
dt 

dM 3 
dt 


= 1 + 2 Ml 


= 1 + 6 M 2 M 3 


dM 4 

dt 


1 + 2Mf + 12 M 2 M 4 


(98) 


Solving 
M 2 = 
M 3 = 


m 4 


these equations one obtains 
tanT 

tanT + | (tanT ) 3 

71 

(cos7) 6 + 5 tanT + 8 (tanT ) 3 + 3(tanT) 
8\/2 


5 


(99) 


where we again used the shorthand notation T = ty/2. 
All the moments diverge at T c = ir/2. This determines 
the gelation time: T c = 7 r/ 2 , or t c = Tt/y/&. Using (99) 


we establish the divergence of the moments 


M 2 1 

M 3 ~ i (t c - t)~ 3 (100) 

, . Zlvy/2 , . 

near the gelation time. Generally 

M n ~ {t c - t )-«(n-i)/2 ( 101 ) 

for the n th moment. Therefore the moments of the cluster 
mass distribution exhibit a multi-scaling behavior. 

The divergence of the moments (101) and multi-scaling 
also arise in the un-driven case dr and they can be un¬ 
derstood using similar arguments. The idea is to use (74) 
with A = 2. The auxiliary time variable is now 


i 

= / dt' M 2 (t 1 ) = — - ln(cosT) 

Jo 2 


( 102 ) 


Since the auxiliary time diverges as T —► T c = tt/ 2, we 
can use continuum approaches; in particular, we can rely 
on the partial-differential equation (74). As explained in 
Ref. the cluster mass distribution attains (for large 
masses) a log-normal form 


»( T ) 


e T/4 m 3 / 2 exp 


(In to)' 


4 r 


(103) 


As a consistency check, let us compute the mass density. 
We write k = In to and find 


E 

m> 1 


IJlCr, 


o-r/i 


dn exp 


21 


K 

4 T 


const 


as r —> oo. The mass density is M = t, so it does ap¬ 
proach a finite value at the gel point: M —► t c = n /y/% 
Similarly, the n th moment grows as 


m>l 


m c m 


0 -t/4 


i(n— 1)t 


da exp 


1 


k \ n — 


21 


K 

4 r 


This reduces to (101) once we notice [see (102)] that 

(104) 


e 2r = 


cosT “ 71 ~ ^ _1 


The moments with n > 1, and even non-integer mo¬ 
ments with n > behave according to (101). The zeroth 


moment, i.e., the total cluster density, exhibits a spe¬ 
cial behavior. This becomes evident once we note that 


J2m> 1 to 3 ^ 2 converges. Therefore N 
we conclude that 


o~t /4 


PH ) 


N(t) ~ (t c - P 8 


Using 


(105) 


Vanishing of the total cluster density at t c shows that 
gelation is complete, namely the entire mass of the sys¬ 
tem is suddenly transformed into gel. 
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V. EXCHANGE PROCESSES DRIVEN BY A 
LOCALIZED SOURCE 

In previous sections we considered mass exchange pro¬ 
cesses supplemented by homogeneous input. Here we in¬ 
vestigate the behaviors when input is spatially localized. 
The densities c m (r, t) then evolve according to a system 
of partial differential equations 


behaviors for the generalized product kernels. At the 
end of this section we show that scaling approaches still 
apply to models in which diffusion coefficients vary with 
mass algebraically (specifically, we discuss models with 
D m ~ m _1 ). 

A. Product Kernel 


da 


— ''y \ CjCj [<^m,z-(-l T 1 2<5 m y] 




T + JS rn ^iS(Y^ 


(106) 


The first term on the right-hand side of Eq. (106) ac¬ 
counts for exchange events. In writing the second term 
we have assumed that the transport mechanism is dif¬ 
fusion; the diffusion coefficients D m generally depend 
on the mass m. The last term represents the local¬ 
ized source, namely the constant flux of monomers (with 
intensity J ) at the origin. It suffices to model flux 
with delta function <5(r) since we are interested in dis¬ 
tances greatly exceeding the size of the region where the 
monomers are injected. 

An infinite system of coupled non-linear partial dif¬ 
ferential equations (106) possesses an outstanding chal¬ 


lenge. Similar equations have been studied in the context 
of aggregation with a localized source [321434] , and mod¬ 
els with mass-independent diffusion coefficients proved to 
be much more amenable to analyses. For exchange pro¬ 
cesses with spatially varying densities, the situation is the 
same. One can already suspect this by noting that when 
D m = D = const, the mass density satisfies a closed par¬ 
tial differential equation, the diffusion equation with a 
localized source: 


dM 

~dt 


DW 2 M + J6(r) 


(107) 


Hence the mass density is known and totally independent 
on the mass exchange: 


M(r,t) = J [ 
Jo 


dt' 


0 (47t Dt') d / 2 


exp 


4 Dt' 


(108) 


in d dimensions. 

We limit ourselves to the physically most relevant 
three-dimensional case. In three dimensions, the mass 
density is asymptotically time-independent: 


M = 


J 

AnDr 


(109) 


More precis ely, t he mass density is described by the 
Coulomb law (1091 when r -C '/Dt. In the following 
we consider physically most important three-dimensional 
case and focus on the region r <C \JDt where not only 
the mass density but all cluster densities are stationary. 

We now analyze two basic reaction kernels, the product 
kernel and the constant kernel, and then mention chief 


For the product kernel K z j = Kij, the governing equa¬ 
tions in the stationary regime read 


0 = ^[(m + l)c m+ i - 2mc m + (to - l)c m _i] 

+ V 2 c m + — <5 m ,i<5(r) (HO) 


where 


JK 

AnD 2 


( 111 ) 


In this section, dimension-full quantities are used (apart 
from the mass of monomers which is still set to unity). 
The densities have dimension of the inverse volume, so 
[ c m ] = L 3 when d = 3. The dimensions of the main 
parameters of the problem are 

[J] = T -1 , [D\ = L 2 T _1 , [K] = L 3 T~ l (112) 


Thus cr has dimension of the inverse length: [<j\ = L — . 

For large masses we can treat m as a continuous vari¬ 
able. The density c = c m (r) satisfies 


a d 2 (mc) 
r dm 2 


V 2 c = 0 


(113) 


Balancing terms gives ^ ~ A, so the typical mass in¬ 
creases linearly with distance: m ~ err. This together 
with the Coulomb law (1091 for the mass density lead to 
the scaling form for the cluster mass distribution: 

TO 


Cm{r)= ID ^{ar) 3 F(/z), 


A 4 = 


(114) 


The scaling form (114) is compatible with the Coulomb 
law (109) when 


/»oo 

/ dmiF(fj.) = 1 
Jo 


(115) 


One can re-write (1141 in a manifestly dimensionless 
form 



D 

Ka 


K)- 3 fV) 


In this and many other equations the distance r from the 
origin naturally appears as a dimensionless combination 
err. The dimensionless parameter 


D A D 3 

- = 47T - 

Ka JK 2 
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provides the proper measure of the strength of the flux. 
By inserting (114) into (113) we obtain 


{li 2 + n)F" + (6/x + 2)F + 6F = 0 
The general solution of this equation is 

C\ + C 2 + 2 In fi — // 


(116) 


Balancing the terms of this equation yields ~ yr, 
leading to m ~ yj'NKr 2 / D. Another estimate of the 
typical mass is m ~ M/N. Balancing these estimates 
and using (109) we get 


er 2 / 3 r _4 / 3 , 

K 


m ■ 


rl/3^1/3 


(124) 


n») = 


(M + l ) 3 


This scaled mass distribution has the small mass singu¬ 
larity unless C 2 = 0. Equation (1151 fixes another ampli¬ 
tude: C\ = 2. Thus F = 2(/x + 1) 3 and (114) becomes 


These estimates determine an appropriate scaling form 
for the cluster mass distribution, e.g., as a scaled mass 
one should use m/(ar) 1 ^ 3 . Putting such a form into (123) 


c m (r ) = 


2Da 2 


and solving the resulting ordinary differential equation 
one finds the scaled mass distribution. The results are 
less cumbersome if we write c m (r) as 


K (to + trr) 3 
In particular, the density of monomers is 

2D 


c,(r) = —r 


-3 


(117) 


(118) 


Cm(r) = (or) 5/3 


H = 


m 


(ar) 1 / 3 


where instead of a defined by Eq. (|111|) we use 
36cr 


9 JK 


a = 


The total cluster density is also given by a neat formula 

D 


N(r) = — r 
v ; K 


(119) 


7r ir 2 D 2 

The scaled mass distribution then reads 

4 /x 


The asymptotic behaviors (118)-(119) are valid far 
away from the origin, more precisely when r ct _1 . 


= 9 {P ? 

The density of monomers is 


1)3 


This is already evident from (117). Another way to see 


this is to consider the governing equation for N (r) which 
is obtained by summing (110) for all to > 1 : 


0 = ——ci + V 2 A r + S(r) 
r D 


( 120 ) 


/ \ 4Z) _ 2 

ClW = 9X r 

and the total cluster density is 

, D a 2 / 3 

^ ~ 9I< r 4 / 3 


(125) 


(126) 


(127) 


(128) 


(129) 


One immediately verifies that 11181) (11191) agree with Near the origin the to tal d ensity is N = J/(4nDr) which 


(120). Near the origin, the source term dominates over again cro sses over to (129) at r ~ a 


reaction, so the total density is N = J/iAitDr). Com¬ 
paring this result with (119) we see that the crossover 
indeed occurs when r ~ a 


Using (129) we find an estimate for the total number 


■=T 


The total cluster density is given by (119) up to dis- 


of clusters in the system: 

D 


-,v 2 > 3 {Dtf ' 6 


tance of the order of y/Dt and it quickly vanishes for 

larger distances. Hence (N(f) ~ f 0 dr Anr 2 N(r) esti¬ 
mates the total number of clusters in the system. Using 
(1191 we obtain 


K 


C. Generalized Product Kernels 


(130) 




B. Constant Kernel 


( 121 ) 


In the continuum limit the density c = c m (r) satisfies 

+ V 2 c = 0 (131) 


M\K d 2 (m x c) : 2 


D dm 2 

The same arguments as in the previous subsections lead 
to the scaling form 


For the constant reaction kernel, J-Qj = K , the gov¬ 
erning equations in the stationary regime become 

0 = NK[cm+i — 2c m + c m _i] +HV 2 c m + JS m! iS(r) (122) 

In the continuum limit the density c = c m (r) satisfies 

NK d 2 c 


Cm(r) = ^ ^ F(n), n= m 


K (ot) 1+2A 


(■arY 


(132) 


where A = 1/(3 — 2A). We omit a rather complicated so¬ 
lution for the scaled mass distribution and only mention 
the asymptotic for the total cluster density 


D dm 2 


V c = 0 


(123) 


D a 


1—A 


N(r) ~ — 


K r 1+A 


(133) 
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and for the total number of clusters in the system 


-£-(<r 2 £>i) l-A/2 
Kg 


The monomer density is given by exactly the same for¬ 
mula (118) as in the model with product kernel and mass- 


(134) 


independe nt dif fusion coefficients, but other densities dif¬ 
fer. Using (138), we compute the mass density 


Equations (132) (H34fo provide a proper description 


only if the total number of clusters diverges in the t —> oo 
limit. This together with (134) give A < 2 implying that 




2 Da 2 
K 


dz 

z 


(139) 


the index A of the generalized product kernels must obey An estimate r ~ \jDt/m max gives m. max ~ Dt/r 2 . Thus 


A < |. Indeed, the rate equation framework which was 
used in the derivations is deterministic, so its applica¬ 
bility is based on the assumption that the total number 
of clusters in the system is infinite. For systems with fi¬ 
nite number of clusters, there are always fluctuations. In 
the exchange processes with a localized source the total 
number of clusters is of course finite, but since we are in¬ 
terested in the large time limit the emergent results are 
consistent only when lim^oo AT(t) = oo. 


D. Models with D m = m 

The assumption that the diffusion coefficients do not 
depend on mass was crucial in the preceding analysis. 
The diffusion coefficients generally decrease with mass, 
however. The decay is usually algebraic, D m ~ m - "; 
see ;35j for arguments in favor of specific values of the 
mobility exponent (e.g. v = 1 and v = 3/2 in some prob¬ 
lems involving two-dimensional clusters). It turns out 
the models with mass-dependent algebraically decaying 
diffusion coefficients are still tractable. The idea is to 
assume the emergence of a stationary mass distribution 
and the validity of scaling. Whenever the results are con¬ 
sistent, one hopes that they are asymptotically exact. 

For concreteness, let us look at a set of models with 
diffusion coefficients inversely proportional to the mass, 
D m = D/m. For such models, the mass density varies as 


^ =D\7 2 N+J6( r) 


(135) 


In the long time limit the left-hand side of Eq. (135) tends 
to zero, so the cluster density is 


N = 


J 


AirDr 


(136) 


(Recall, that we consider the three-dimensional case.) 
We know N(r), so the model with constant kernel is the 
simplest. In that model, the densities satisfy 

0 = — [c m +i — 2c m + c m - 1] + V--—h — S m i<5(r) (137) 
r mD 

The scaling solution is found using the same techniques 
as previously to yield 


c m (r) = 


m 


2Da 2 


K (to + err) 3 


(138) 


M(r,t) = 


2 Da 2 
K 


In (a 2 Dt) — 31n(or) — - 


(140) 


Models with D m = D /to and non-constant kernels are 
more complicated. Let’s try, however, to probe chief be¬ 
haviors relying on scaling and assuming the emergence of 
a stationary mass distribution. For models with general¬ 
ized product kernel we arrive at a familiar scaling form 


D a 2 

c ™( r ) = yy 


K (err) 


II = 


(or)' 5 


(141) 


with S = 1/(1 — 2A). Even the the mass distribution is 
stationary, the mass density may be non-stationary as it 
is already exemplified by Eq. (139) in the case of A = 0. 


Using the scaling form (141) we arrive at an estimate of 
the mass density 


„ ,, , D 2 — 2X 

M{r) ~ —u 1 - 2X 
K 


(142) 


The dependence of M(r) on the distance suggests that 
the mass density is stationary for A < 0. Using mass 
conservation, JJ” 1 " 1 dr r 2 M (r) ~ Jt, we conclude that 
(142) is valid for r < r max with 


err u 


(a 2 Dt) 1 


(143) 


and the mass density quickly vani shes for r > r max . 

A more rigorous derivation of (142) in the A < 0 re¬ 


gion could be obtained if we determine the scal ed mass 
distribution. By inserting the scaling form (141) into 


M\K d 2 (m x c) 


D 


dm 2 


-X7 Z — = 0 

TO 


(144) 


we obtain 


^ liF" + (2 - A )F' + 1 - A F + Ix(li X F )" = 0 (145) 

2 (i 

with I\ = |(1 — 2A) 2 / 0 °° d/i ii x F(il). Apart from the al¬ 
ready solved case of A = 0, Eq. ( |145[ ) apparently does not 
have solutions in terms of known special functions. Still, 
the large mass limit can be probed analytically. In this 
limit, the last term in Eq. (1451 is negligible. Dropping 


this term and solving the resulting equation we obtain 

^- 2 (i-A) for ( 146 ) 


F 


For A < 0, the integral f°° dji liF(fi) conv erges and hence 
the mass density does scale according to (142). 
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When 0 < A < the mass density is non-stationary: 

/* Mmax 

M{r,t) ~ ^(crr)^ 1 

(H 5 - 1 


2 » 

Da 2 
K 

Da 2 4 
A' 

Her 2 /Dfo 2A 


dp pF(p) 

3 

f £)t/[r 2 (or)' 5 ] 


dp /i 


2A—1 


A' 


(147) 


where we have used a n estimate TO max ~ Dt/r 2 . Also for 
0 < A < \, Eq. (1471 is valid for r < r max with 


OTn 


{a 2 Dt) 


(1 —2A)/(3—4A) 


(148) 


Us ing (136) together with the estimates (143) and 
(1481 of the maximal distance we find that the total num- 

N ~ {a 2 Dtf (149) 


ber of clusters scales as 


Ka 


with 


/3 = 2(1 — 2A) x 


(3-2A)- 1 A < 0 
(3-4A)" 1 0 < A < 


The above predictions are inconsistent when A > i. 
The most natural guess is that the cluster mass distri¬ 
bution does not reach a stationary regime in this region. 
It would be interesting to develop a quantitative under¬ 
standing of the behavior in the A > | region. 


VI. CONCLUSIONS AND DISCUSSION 

We studied mass exchange processes driven by a source 
of monomers. We considered both homogeneous and lo¬ 
calized input. For processes with homogeneous source, 
we employed scaling methods and analyzed in detail 
three basic exchange kernels—the product kernel, the 
sum kernel, and the constant kernel; we also probed a 
one-parameter family of generalized product kernels. We 
also exactly solved the mass exchange process with prod¬ 
uct reaction kernel supplemented by the homogeneous 
source of monomers. In all these systems, a steady state 


is never reached. A generalization of our approach to 
exchange processes with kernels K l:j = i a j b + i b j a and 
homogeneous input seemingly leads to the same conclu¬ 
sion. It would be interesting to find a general proof that 
steady states are impossible in exchange processes driven 
by a homogeneous source. 


We showed that the mass exchange processes with spa¬ 
tially localized input can be simpler than the correspond¬ 
ing processes with homogeneous input. This is obvious 
e.g. by comparing the cluster mass distributions for the 
exchange process with product kern el: E qs. (25) and (31) 
for the homogeneous input and Eq. (117) for the localized 
input. The chief reason is that the cluster mass distri¬ 
bution is often stationary when the source is spatially 
localized. We demonstrated that stationary mass dis¬ 
tributions naturally emerge for exchange processes with 
mass-independent diffusion coefficients, e.g., for exchange 
processes with generalized product kernel, K, 3 = ( ij) x 
with A < |, where the bound on the homogeneity index 
assures that the deterministic rate equation description 
is asymptotically relevant. 


We outlined possible behaviors for the mass exchange 
processes with spatially localized input in the situations 
when diffusion coefficients depend on mass. Specifically, 
we analyzed models in which diffusion coefficients are in¬ 
versely proportional to the mass. We presented evidence 
that the cluster size distribution may not reach a station¬ 
ary state. This depends on the details of the exchange 
kernel; e.g., for exchange processes with generalized prod¬ 
uct kernel a stationary state is apparently never reached 
when A > |. The behaviors in the continuously evolving 
regime are the challenge for future work. 


The power of scaling in the context of exchange pro¬ 
cesses is impressive. Another general lesson is that seem¬ 
ingly intractable systems of infinitely many coupled non¬ 
linear partial differential equations for spatially varying 
densities can be simpler than their homogeneous counter¬ 
parts. This is not clear a priori, and even the emergence 
of the stationary distribution does not make the spatial 
problem manifestly simpler as it involves the Laplace op¬ 
erator rather than the first derivative in time. 


Finally we mention that exchange processes often oc¬ 
cur in low spatial dimensions where the rate equation 
approach is invalid. It would be interesting to analyze 
simplest exchange processes in one dimension. 
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